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Abstract. The main objective of this article is to postulate a principle of 
interaction dynamics (PID) and to derive unified field equations coupling all 
four forces. PID is a least action principle subject to divyi-free constraints 
for the variational element with A being gauge fields. The unified field equa- 
tions of the coupled interactions of four forces are then uniquely derived based 
on 1) the Einstein principle of general relativity (or Lorentz invariance) and 
the principle of equivalence, 2) the principle of gauge invariance, and 3) the 
PID. The unified model gives rise to a complete new mechanism for sponta- 
neously breaking gauge-symmetries and for energy and mass generation, which 
provides similar outcomes as the Higgs mechanism. For the electromagnetic 
and weak interactions alone, we derive a totally different electroweak theory, 
producing the three vector bosons and Z, as well as three Higgs bosons- 
two charged and one neutral. One important outcome of the unified field 
equations is a natural duality between the interacting fields [g, AjW"" , B^), 
corresponding to graviton, photon, vector bosons and Z and gluons, and 
the adjoint fields , rp^ , , (j>''g) , which are bosonic fields. The interaction 
of the bosonic particle field 3> and graviton leads to a unified theory of dark 
matter and dark energy and explains the acceleration of expanding universe. 
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1. Introduction 

As we know there are four forces/interactions in nature: the electromagnetic 
force, the strong force, the weak force and the gravitational force. There are suc- 
cessful theories to describe these interactions, including the Einstein general theory 
of relativity for gravitation, the quantum electromagnetic dynamics (QED) for elec- 
tromagnetism, the Weinberg-Salam electroweak theory unifying the weak and elec- 
tromagnetic interactions [l]|7l[6], and the quantum chromodynamics (QCD) for the 
strong interactions. In particular, the Standard Model, a SU{3) (E) SU{2) (E) U{1) 
gauge theory, provides an attempt for all known interactions except gravity; see 
among many others [4 . 

This article is an attempt to derive unified field equations coupling all four 
interactions. Apparently the unified field equations will be built upon the success 
of gauge theory and the Einstein general theory of relativity. 

Motivated by the Higgs mechanism of electroweak theory and the new gravita- 
tional field equations incorporating dark energy and dark matter [4j [5] , the starting 
point of the study is to postulate a new least action principle for actions subject to 
energy-momentum conservation constraints, which we call the Principle of Interac- 
tion Dynamics (PID); see Section 3 for the precise statement of this principle. 

The main motivation of postulating PID is as follows. Due to the presence of 
dark energy and dark matter, the energy-momentum tensor Tij of normal matter 
is no longer conserved: 

Therefore by an orthogonal decomposition theorem [5] , the energy- momentum ten- 
sor Tij can be decomposed naturally as 

T^j = T,j - -^-^D.Djip, 
D% = 0, 

where ip is a scalar function defined on the space-time manifold. On the other hand, 
the Euler-Lagrangian of the scalar curvature part of the Einstein-Hilbert functional 
is divergence-free due the the Bianchi identity: 

D\R^J - lg^JR) - 0. 

Hence the modified gravitational field equations become 

1 SttG- 

which are equivalent to 

1 SttG 
(1.1) - -g,jR = —T^ii - AI>i<^- 

The derivation here is equivalent to studying the least action of the Einstein-Hilbert 
functional under divergence-free constraint for the variational element. To see this, 
consider the Einstein-Hilbert functional Leh given by 

Teh — I 
J M 
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Then PID amounts to saying that for any X — {Xij} with D'gXij ~ 0, 

Um ^[LsHigij + AXij) - LsHigij)] = {SLEH{gij), X) = 0. 

By the orthogonal decomposition again, 

SLEnigi]) = -DiDjip, 

which is exactly the gravitational field equations taking into consideration of 

the presence of dark matter and dark energy. 

In order to take into account all four interactions, the natural constraint is then 
the divA-free constraint, leading to the PID. Here A represent the gauge fields. 

The unified field equations of the coupled interactions of four forces are then 
uniquely derived based on the following principles: 

• the Einstein principle of general relativity (or Lorentz invariance) and the 
principle of equivalence, which amount to saying that the space-time is a 
4-dimensional Riemannian manifold {M,gij), 

• the principle of gauge invariance |5 , 

• the principle of interaction dynamics. 

The unified field equations are given by ()4.4p - ()4.10|) . Here are a few main ingredients 
of these field equations: 

First, the Lagrangian action functional L is naturally the sum of three im- 
portant parts as given in (j4.2|) : the classical Einstein-Hilbert functional Ceh, the 
Glashow-Weinberg-Salam functional Cgws without Higgs field terms, and the stan- 
dard QCD functional Cqcd- This Lagrangian action contains no tuning terms and 
is uniquely determined. 

Second, thanks to the div^i-free constraint, the Euler-Lagrangian of the ac- 
tion functional is balanced by gradient fields, resulting the unified field equations 
coupling all four forces. 

Third, the gradient fields break spontaneously the gauge-symmetries. This gives 
rise to a complete new mechanism for spontaneously breaking gauge-symmetries 
and for energy and mass generation, which provides similar outcomes as the Higgs 
mechanism. 

Fourth, one natural outcome of the new mechanism is the adjoint fields 
(fp, (j)^, (/)^. Amazingly, there is a natural duality between the interacting fields 
{g, A, W, B^) and the adjoint fields ($^, (jP , (j)'^ , (j)]^): 



(1.2) 











A, 








w; ^ 




for a ~ 


1,2,3, 






for h ~ 
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On the left hand side of this duality, the interacting fields represent the following 
particles: 

{g^iy} represents the graviton g, 

represents the photon 7, 

W° represents the vector bosons and Z, 

B^^ represents the gluons gb- 

On the right hand side of the duality (|1.2p , each field is a bosonic field: $ represents 
a vector bosonic field, and 0", i/)^, (ffi^ represent scalar bosonic fields. 

Fifth, the combination of the three Higgs bosonic fields on the right hand side 
of the duality induces three Higgs boson particles given by 

with (fy^ having ± electric charges, and with being neutral. Note that the 
classical Weinberg-Salam theory induces one neutral Higgs boson particle. All 
three Higgs bosons deduced here possess masses, generated by the new mechanism. 

In fact, consider only the electromagnetic and weak interactions and ignore the 
effect of other interactions, we derive a totally different electroweak theory. Again 
this electroweak theory produces the three vector bosons and as well as the 
three Higgs bosons. The spontaneous gauge-symmetry breaking is achieved by the 
constraint action without Higgs terms in the action functional. 

Sixth, the new vector particle field $ corresponding to the gravitational field 
{gij} is massless with spin s — 1. This particle field corresponds to the scalar 
potential field, caused by the non-uniform distribution of matter in the universe as 
explained in [5]. In other words, the theory here demonstrates that the interaction 
between this particle field $ and the graviton leads to a unified theory of dark 
matter and dark energy and explains the acceleration of expanding universe. 

Seventh, the field 0", adjoint to the electromagnetic potential A^, is a massless 
field with spin s = 0, and it needs to be confirmed experimentally. We think it is 
this particle field (jp that causes the vacuum fluctuation or zero-point energy. 

Eighth, we conjecture that the combination of the eight scalar bosons 0^ 
corresponding to gluons 56 (1 < < 8) leads to the Eightfold Way mesons: 7r+, tt", 
A'-,i^o, i^",r;0. 

2. Orthogonal Decomposition and Variations with div^-FREE 

Constraint 

In this section, we address 1) the orthogonal decomposition of tensor fields with 
respect to gauge-connections A, and 2) the variations of functionals with div^i- 
free constraint. They are crucial for the derivation of the unified field equations of 
coupled interactions of all forces in the subsequent sections. 

2.1. Orthogonal Decomposition Theorems. In 0, we obtained some orthog- 
onal decomposition results for general tensor fields on a Riemannian manifold. For 
convenience, we briefly recall them in the following. 
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(2.2) 

Equality 

Remark 2.1 

(2.3) 



Theorem 2.1 (Ma and Wang [5 ). Let {AI,gij) be a Riemannian manifold (in- 
cluding the Minkowski type manifolds) without boundary, and u G L'^iT^M) with 
r + s > 1, where L'^{T^'M) is the space consisting of all (r, s)-tensor fields on M 
which are square integrable. Then the tensor field u can be orthogonally decomposed 
into 

(2.1) u = Dip + v, div u = 0, 

for some ip G H^iT^-Hl) or (p e H'^{Tl_^M), i.e. 

I Dp-Vy/gdx^O, g^ det {gij). 

is ensured by the formula 

/ Dip ■ Vy^dx = ip ■ div v^dx. 

Jm Jm 

Theorem 2.1 implies that if a tensor field u G L'^iT^M) satisfies 
u-Xy/gdx = yX e L^{T^M) with div X = 0, 

then there exist a tensor field (p such that 

(2.4) u = Dip. 

Based on the orthogonal properties (I2.ip - (|2.4p . we can derive the follow scalar 
potential theorem. 

Theorem 2.2 (Ma and Wang fE'). A ssume that the first Betti number of M is zero, 
i.e. any loop in M can shrink to a point. Let F be a functional of the Riemannian 
matric {gij}- If {dij} an extremum point of F with the divergence-free constraint, 
i.e. 

5F{g,j) ■ X^dx = ^F{g,, + \X,,)\x=o = 

for all X ~ {Xij} G ^^(Tj Af) and d\VgX — 0, where divg is under the Levi-Civita 
connection of {gij}, then there is a scalar function ip G H'^{M) such that the metric 
{gij} satisfies that 

(2.5) SF{g,j) = D^cp. 

We now generalize the above decomposition theorems to gauge connections, 
which are crucial to establish the dynamical theory of interaction fields. Hereafter 
we still assume that {M, gij ) is an n-dimensional Riemannian manifold without 
boundary, where the metric {gij} includes the Minkowski type Riemann metric, 
i.e., with a proper coordinate system, {gij} can be expressed as in the form 



(2.6) 



-1 





G 



where G is a positive symmetric (n — 1) x (n — 1) matrix. 
Let u G L^iTlM). We define the operators Da and div^ as 



(2.7) 
(2.8) 



Dau = Du + u® A, 
div^M — divM — A - u, 
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where A is a vector or covector field, D and div are the usual gradient and divergent 
operators. It is clear that for Da and divyi we have the integral formula 

(2.9) / D^u ■ Vy/gdx — — u ■ Ayv AV\/gdx. 
Jm Jm 

Then, for the differential operators given by (|2.7p and (|2.8p we have the following 
orthogonal decomposition theorem. 

Theorem 2.3. Let A is a given vector field or covector field, and u G L'^{T^M) 
with r + s > 1. Then u can he orthogonally decomposed into 

(2.10) u — Da^P + v with div^f = 0, 

where ip G H^{T^-Hl) or ip e H^{T^^^M). 

The proof of Theorem [23] is the same as that of Theorem 12.11 see Ma and Wang 
[5], and we omit the details. 

A few remarks are now in order. 

First, Theorem 12.31 is a generalized version of Theorem 12.11 and is reduced to 
TheoremOif ^ = 0. 

Second, when M is a compact manifold and (gij) is a positive matrix different 
from that given in p.6p . then u can be further orthogonally decomposed into 

u = Dap + V + h, 

divAV = 0, DAh — 0, div^/i = 0, 

where h is a harmonic field, and the space 

n:^{he L^{TIM)\ DAh = 0, div^/i = 0} 

is of finite dimensional. 

Third, the orthogonal decomposition in Theorem 12.31 is very important for the 
interaction field theory developed in this article, because the vector fields A in 
(|2.7p and (|2.2I) represent gauge fields in the interaction field equations, and lead 
to a new energy and mass generation mechanism, different from the famous Higgs 
mechanism. 

2.2. Variations with div^-free constraint. According to Theorem 12. 31 we gen- 
eralize the scalar potential theorem (Theorem l2.2[) to the case where the variations 
of functionals F are under the div^-free constraint. 

Let F — F{u) be a functional of tensor fields u. We say that an extremum point 
Uq of F{u) is under the div^i-free constraint, if 
(2.11) 

d f 
{5F{uo), X) = -77-^(^0 + \X)= / 5F{uo) ■ Xy/^dx = Vdiv^X = 0. 

In particular if is a functional of Riemannian metric , i.e. uo = g^j a Riemann- 
ian metric, then the differential operator Da in divy^X in (j2.11l) is 

D\^D' + A' or {DAh^D, + A. 

and the connection in D is respect to the extremum point uq — (?° . 

Then we have the following theorems for the div^-free constraint variations. 
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Theorem 2.4. Let F = F{gij) he a functional of the Riemannian metric. Then 
there is a vector field $ G L^(TM) such that the extremum points {gij} of F with 
the divyi-/ree constraint for vector field A satisfy the equations 

(2.12) SF{g,j) = D<i> + A(g)'i>. 

In particular, if the vector field A — 0, then the conclusions of Theorem \2.B hold 
true. 

Theorem 2.5. Let F = F(u) be a functional of a vector field u. Then there is a 
scalar function G H^(M) such that for given vector field A, the extremum points 
u of F with the dwA-free constraint satisfy the equation 

(2.13) 5F{u) ^ Dip + LpA. 

Proof of Theorems \2.4\ and \2.5[ First, we verify Theorem 12.41 By (|2.11|) . the ex- 
tremum points {gij} of F with divA-free constraint satisfy 

/ SF{g,j) ■ Xy/^dx ^0 yX e L^{T^M) with divAX ^0. 
Jm 

It imphes that 

(2.14) 5F{gij)LL]j{T^M) = {w G L^{T^M) \ div^u = 0}. 

By Theorem 12.31 L^(T2 M) can be orthogonally decomposed into two parts 

L^{T^M) = Ll){T^M) ® G^{T^M), 
G^{T^M) = {DaB\ B e H\T^M)}. 

Hence, it follows from (|2.14p that 

6F{g,,) e G^{T^M), 

which means the equality (I2.12p holds true. 

We are now in position to prove Theorem l2.5l Similarly, for the extremum vector 
fields u oi F with div^-free constaint satisfy 

(2.15) [ SF{u) ■ Xyf^dx = yXeL^{TM), div^X 
Jm 



- 0. 

Im 

In addition. Theorem 12.31 means that 

L^{TM) = L%{TM) ® G^{TM), 
L^{TM) = {ve L^{TM)\ divAV = 0}, 
G^{TM) ^ {DAf \ V e H^{M)}. 

We infer from ([235]) that 

5F{u) e G^{TM). 

Theorem 12.51 follows, and the proof is complete. □ 
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3. Principle of Interaction Dynamics 

We know that there are four fundamental forces in nature: the gravitational 
force, the electromagnetic force, the weak and the strong forces. There are success- 
ful theories to describe these interactions, including the Einstein general theory of 
relativity for gravitation, the QED for electromagnetism, the Weinberg-Salam elec- 
troweak theory, and QCD for the strong interactions. In particular, the Standard 
Model, a SU{3) (X) SU{2) ® C/(l) gauge theory, can rather successfully describe aU 
known interactions except gravity. 

However, there is still no ideal unified theory to explain all four forces. For 
this purpose, we propose a basic principle, which we call Principle of Interaction 
Dynamics (PID). This principle is based on the variational theory developed in the 
previous section, and can be considered as a revised version of the Lagrangian least 
action principle for interacting fields. 

Principle of Interaction Dynamics (PID). For all physical interactions there 
are Lagrangian actions 

(3.1) L(g,A,V)= / L{g,j,A,i,)^dx, 

Jm 

where g — {gij} is the Riemann metric representing the gravitational potential, 
A is a set of vector fields representing the gauge potentials, and ^ are the wave 
functions of particles. The actions iS. 1\) satisfy the invariance of general relativity 
(or Lorentz invariance) and the gauge invariance. Moreover, the states {g,A,^p) 
are the extremum points of iS.l]) with the divyi-/ree constraint. 

According the above principle and Theorems 12.41 and 12.51 the field equations 
with respect to the action (|3.ip are given in the form 

(3.2) — -L(5,A,V) = (A + ^afe4')*.. 

^ fe=i 

f) ^ 

(3.3) jj^L{g, A, V) = (A + P^^oA't)^'^ 

(3.4) _L(g,A,^)=0, 

where A^, - ■ ■ ,A^ are the gauge vector fields for the electromagnetic, weak, and 
strong interactions, $ = ($i,... is a vector field induced by gravitational 

interaction, ip'' are scalar fields generated from the gauge field A'^ = (Aj^, . . . , A^), 
and ckfe, Pkj (1 ^ k,j < N) are coupling parameters, which may depend on physical 
quantities as energy densities. 

Now, we need to give some explanations for these field equations p.2p - p.4p . 

1 . As explained in the Introduction, there are several reasons behind the div^- 
free constraint. First, the coupling energies of interactions do not usually 
possess variational structure. In other words, the terms in the right-hand 
sides of (|3.2p and (|3.3|) cannot be put into the Lagrangian. Second, the con- 
straints eliminate the extra freedoms induced by the symmetry-invariance 
of Lorentz or general relativity and the gauge invariance. Hence each of (p'^ 
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determines a gauge. Third, the div^-constraint can be regarded as impos- 
ing conservations of energy-momentum. 

2. The equations p.2p - (|3.4l) preserve covariance of Lorentz or general relativ- 
ity. However they break the gauge symmetry, caused by the constrained 
variational principle. In fact, it is utterly important that the field equations 
must break the gauge symmetry to avoid the existence of infinite number 
of solutions. However the action functional must preserve the gauge sym- 
metry. Hence the least action must be taken with constraints, leading 
naturally to spontaneous gauge-symmetry breaking. 

3. The spontaneous gauge-symmetry breaking caused by the terms afeAf $j 
and f^ijA'^if'^ in (|3.2I) and (|3.3I) leads to energy creation. Also, for the 
weak interaction, the gauge-symmetry breaking gives rise to a new mass 
generation mechanism, different from the classical Higgs mechanism. 

4. The fields g, A, ip together with $ and tp'^ as unknowns of p.2p - (l3.4p char- 
acterize the interacting states. The fields $ and (p'^ are called the adjoint 
fields of g and A'^ respectively. 

5. The coupling parameters ak represents the coupling strength between the 
gravitation and the gauge field A'^ , and /3k j represents the coupling strength 
of A'' and A^ . They are determined by physical experiments and physical 
laws. These parameters are natural and required due to the energy-coupling 
of different fields. 



4. Unified Field Equations Coupling Four Forces 

In this section we derive a unified model for all four interacting forces, based 
on the Principle of Interacting Dynamics. The action functional is the natural 
combination of the Einstein-Hilbert functional, the electroweak action, and the 
QCD action, and is given in the following form; 



(4.1) 




where 




(4.2) 



6 



+ X! {ifc(«7^^M - mk)Lk 



+ 



^QCD 



ml)qk, 



where R is the scalar curvature of spacetime Riemannian manifold, S is the mo- 
mentum density, W° (a = 1,2,3) are the electroweak gauge potentials, is the 
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electromagnetic potential, B'^^ (1 < 6 < 8) are the color gauge potentials, and 

= D,Bl - D^Bl + g2g''''BlBt 
D^Lk - (V^ - i^W^aa + 

Df,qk = (V^ + i^Bl\b)qk, 

where gi and g2 are constants, e is the charge of an electron, f'^'"^ and g''^'' are 
the structure constants of SU{2) and SU{i) respectively, aa (a = 1,2,3) are the 
Pauli matrices. A;, (1 < 6 < 8) are the Gell-Mann matrices, Lk = {tpi^k i V'fc )' are the 
wave functions of left-hand lepton and quark pairs (each has 3 generations), are 
the right-hand leptons and quarks, qk — {qki,qk2,qkzY is the k-th flavored quark. 
Dp is the general covariant derivative, and is the Lorentz Vierbein covariant 
derivative [1]. 

For the gauge fields V^, we have 



By (|3.2p and (|3.3p . it follows from (|4.3p that the field equations coupling all four 
interactions are given as follows: 

1 SttG 

(4.4) R^..~^g^.R+—T^,. 



(4.5) d^'{df,A,-d,A^)-eJ2 



{D^ + aoAf, + awaWl + aBjBl)<^^ 
1 



6 r , 



k=l 

= {D, + l3oA, + PwaW^ + fiBjBi)4P, 

3 6 

(4.6) - I ^ f^'-g'^'^Wl^.W^^ + I ^ Lk^.cjaLk 

6=1 fe=l 
A;=l fc=l 

= (D. + + Cl^^ + O'biBD^'b, 

(4.8) (i7^5p - mfe)Lfc forfeit, 2, 3, 

(4.9) {i^D^ - mi)V^ = for fc = 1, 2, 3, 

(4.10) (z7^5p - ml)qk. = for A: = 1, • • • , 8, 
where 

6S c"^ 
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Of course, the above field equations need to be supplemented with coupled gauge 
equations to fix the gauge. We think the gauge equations should be given by 

(4.11) Df'Af, = constant, D^W^ = constant, D^'B^ = constant. 

From the field equations (|4.4p - (l4.7p . we see that there is a corresponding relation 
between the interacting fields {g, A, W"", B^) and the adjoint fields ($^, i/)^, t/)^) 
given as follows: 



(4.12) 











A, 










^ rw 


for a ~ 


1,2,3 




^ 4>'b 


for h = 





In addition, we know that 

{9fii^} represents the graviton g, 
A^i represents the photon 7, 
W° represents the vector bosons and Z, 
B^ represents the gluons gb- 
Taking divergence on both sides of (|4.4p - (|4.6p respectively, we obtain that 
D^D^* + (divyO* = 0(1), 

+ (divF2)./)° = o(l), 
D^'D^cj)^ + {divVs^)rw = 0(1) for a - 1, 2, 3, 

Df^D^^'s + (divF3'')</)^ - 0(1) for & = 1, • • • , 8, 

which imply that $ represents a vector boson, and (jP, 't'^B represent scalar 
bosons. Here 

Vi = ao^M + awaWj^ + asjBf^, 
V2^^oA,+l3waW^ + l3B,Bi, 

V3 — I^qAu + K^fjVF^ + KgjBl, 

vl = eiA, + e\^,wt + e'^^Bi. 

In other words, each of the fields {g^ny}, W° (a = 1,2,3) and i?^ (6 = 
1, • • • ,8) corresponds to a bosonic field. Therefore the corresponding relation (|4.12p 
leads to the following conjectures: 

(1) There are three Higgs boson particles given by 

with (fy^ having ± electric charges, and with 0^ being neutral. Note that 
the classical Weinberg-Salam theory induces one neutral Higgs boson par- 
ticle. All three Higgs bosons deduced here possess masses, generated by 
the new mechanism here, different from the classical Higgs mechanism. 

(2) The combination of the eight scalar bosons 0^ corresponding to gluons gb 
(1 < 6 < 8) leads to the Eightfold Way mesons: 7r+,7r~, tt°,K+, K-,K°, 
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(3) The new particle field $ corresponding to the gravitational field {gij} is 
massless with spin s = 1. This particle field corresponds to the scalar 
potential field, caused by the non-uniform distribution of matter in the 
universe as explained in |5j. The interaction between this particle field $ 
and the graviton leads to a unified theory of dark matter and dark energy 
and explains the acceleration of expanding universe. See also Conclusion l5.ll 
below. 

(4) The field cfP, adjoint to the electromagnetic potential A^, is a massless field 
with spin s = 0. Of course, the particle adjoint to photon 7, still needs 
to be confirmed experimentally. We think it is this particle field (p'^ that 
causes the vacuum fluctuation or zero-point energy. 

(5) From the mathematical point of view, the classical Einstein field equations 
contain 10 gravitational field equations solving for 6 unknown functions in 
gij due to coordinate transformations, leading to an ill-posed problem. This 
ill-posedness is resolved in the field equations (|4.4p - (|4.10p due to the pres- 
ence of the vector potential derived from the variation with constraint. 



5. Gravitational Field Equations 

We know that the action of gravitation is the well known Einstein-Hilbert func- 
tional given by 

(5.1) L= [ 

J M 

When the other three interacting forces are very weak, the gravitational field 
equations decouple with the gauge fields A = {A^,--- In this case the 

adjoint field $ of the metric g in (|3.2p becomes a gradient field of a scalar function 
as demonstrated by the authors in [5]: 

$ = -Dip. 

Namely the decoupled gravitational field equations for (|5.ip are taken in the fol- 
lowing form 

1 SttG" 

(5.2) i?y - -g,jR = —^ii - DiDjtp. 

The interested readers are referred to [5 for more detailed discussions on the above 
modified gravitational field equations and their physical implications. In particular, 
the new field equations give rise to a unified theory for dark energy and dark matter, 
and provide an explanation for the Rubin rotational curves and for the acceleration 
of expanding galaxies. 

In the case where the gauge fields A can not be ignored, the field equations p.2p 
are expressed as 

1 SttG ^ 

(5.3) - -g,,R = -—^T^J + + ^ XkA^^j 

k=l 

Since the strong and weak interactions are short range forces, in the large scale 
domains, only the electromagnetic field A^ = {Aq, Ax, A2, A^) is retained in (|5.3p . 
Hence the large scale coupling interacting gravitational field equations are of the 
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following form 

1 SttG 

(5.4) R^^ - -g^^R = — ^^A"' + ^t^'^'' + •^^p*!': 
where A is a coupling parameter, and 

(5.5) $^ = ($o,«'l,*2,«'3) 

is the adjoint vector potential of gravitation, which is regarded as the effects induced 
by the coupling of matter and electromagnetic fields. 
A direct consequence of (I5.4p - (|5.5p is as follows: 

Conclusion 5.1. There exists afield different from photon field, which is a vector 
boson $^ as in 115. 5\) . possessing zero mass and spin s — 1. 

To see this, let the momentum tensor T^i^ = (in vaccuum), and taking diver- 
gence on both sides of (|5.4p . by div(i?^i, — ig^^i?) = 0, we have 

(5.6) D'^Df.'i',, + A(D^yl^$^ + ■ L^^^^) = 0. 

Since the fields and are weak and the electromagnetic gauge df^A^ — , we 
can ignore the second order terms and infer from equations ()5.6p that 

(5.7) L>^L>^$ = 0. 

In weak gravitation, the operator D^Df^ is approximately the Klein-Gorden oper- 
ator. Hence (15.71) implies Conclusion 15. II 

As concluded in the previous section, the massless with spin s = 1 particle field 
$ corresponds to the scalar potential field, caused by the non-uniform distribution 
of matter in the universe as explained in [5]. The interaction between this particle 
field $ and the graviton explains the dark energy and dark matter. 

6. Electroweak Theory 

6.1. Weinberg-Salam Model. We first recall the classical Weinberg-Salam elec- 
troweak theory, which is a SU{2) x C/(l) gauge theory. The action density is given 

by 

(6.1) C^Cg + Cf + Ch, 

where Cq is the gauge part, Cf the fermionic part, and Ch is the scalar Higgs 
sector. They are given as follows [31 [2]: 

(6-2) Lf = iR-z^'D^R + iL-f^'D^L, 

Ch = D^,4)'^D^'4) - A(0V - af + Ge{L(t)R -f i?0^L), 
where 7^ = (7^,7^,7^,7'^) are the Dirac matrices, and 

w% = d,w: - d,w; + gif'^'miwi, 

Ffil, = dfj_B^ — d^B^, 
D^,R = [d^ + ig2B^)R, 

D^L^{^^ + ^^B^-^^W^^<7k)L, 
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where gi and 32 are the couphng constants, Z*^'^ (1 < k,i,j < 3) are the structure 
constants oi SU{2), ak (1 < k < 3) are the Pauh matrices, and 

<■ = {W^M.W^M) for fc = 1, 2, 3, 
= {Bo, Bi, B2, B3), 



L=\ I , i? = V^, 



R 

e ' 




where W!^ is the Yang-Mills gauge field corresponding to the k-th generator of 
SU (2), Bf^ is the gauge field with respect to f7(l), ■i/'i' and tp(^^ are the wave functions 
of left-hand electron and neutrino respectively, tp^ is the wave function of right- 
hand electron, and 4> is the Higgs scalar field. 

The complexity of the Weinberg-Salam model (16. ip and (16. 2p is due to the cou- 
pling to the Higgs field (f), which breaks the gauge symmetry, and creates mass. The 
basic idea, called the Higgs mechanism, is as follows. The action ()6.ip is invariant 
under the SU{2) gauge transformation: 

g-(i/2)e.-<T,^^ 
(6.3) R^R, 

gi 

Bfi ^ Bf^, 
and the U{\) gauge transformation: 

L e^'^^"^L P a real function, 

R e'^R, 

(6.4) 

B^^B^ + —d^/3. 
92 

Hence the variational equations 

(6.5) ^L{W,B,L,R,cj>)=0 
are covariant for the action 

(6.6) ^ " / i^G + CF]dx. 
However, (|6.5p are of the following form 

(6.7) d^Wj:^ = o{W) the higher order terms of W, 

which implies that the field particles described by W'' are massless. This is a 
contradiction to the physical fact that the W'^ particles have mass. 
If one adds a term 
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into the action Cg-, then the action 

breaks the symmetry for the gauge transformation (j6.3p and (j6.4l) . 

Higgs mechanism provides a solution when we add the Higgs action Ch into 
(fO]) . We infer then from dSZ]) and that 


a 

is an extremum point of (|6.ip . Hence 

$o = (O,O,O,O,0o) 

is a solution of 

SL = 0, L = J[Cg + Cf + CH]dx. 
Consider the translation for $ = (VF, B, L, R, (p): 

$ = $ + $0, $ = (W?, 5,1,^,0). 
Then the variational equations of L for $ are given by 

(6.8) = ^ \ +oiW,B)=0, 

where M is the mass matrix. The equations (|6.8p show that the particles {W„, B,y) 
have masses due to the spontaneous gauge-symmetry breaking for the gauge trans- 
formations (|6.3I) and (|6.4I) . 

In the Weinberg-Salam model, there are three mass bosons Z^, W^, and 
photon Afj_ induced from VF^ and i?^ as follows: 

V 9i + 92 

vat + gi 



where the Weinberg angle Ow is defined by 



a 31 -a 92 



The values of the Weinberg angle 9w and masses Mz , My/ are determined as 



sin^ Ow = 0.2325 ±0.008, 



(6.9) Mz = 91.173 ± 0.02 GeV/c^, 

Mw = 80.22 ± 0.26 GeV/c^. 
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6.2. Electroweak theory based on Principle of Interaction Dynamics. 

There are two shortcomings for the classical Weinberg-Salani model. The first is 
its complexity, and the second is that the electromagnetism can not be decoupled 
in a natural and apparent manner. 

In this section, we present a much simpler electroweak theory. 

We take the action 



(6.10) 



where 



1 



1 



+ L{iJ^'D' - m)L + Riijf'Df, - me)R 



da Ay ~ di/Aii 



D'„ = d„, + ieA„, + \Wl:ak, 



ieA,,. 



Obviously, for the case where there is no weak interaction present, VK^^ — — 0, 

and the action (|6.10p reduces to the classical QED action: 



L. 



QED 



dx. 



By (|3.3p and p.4p . from (|6.10p we obtain the electroweak interacting field equa- 
tions as follows: 



(6.11) 

(6.12) 
(6.13) 
(6.14) 

where 



1=1 

= {d, + akWl; + aoA,)(j)'l for fc = 1, 2, 3, 
d^Ff,, + 4 + = (5. + PjWi + /3o^.)</'2, 
{il^D^, - m)L = 0, 
{ij^D^ - A'QR = 0, 



4 = iXLYcTkL, 
Jl — ieLj'^L, 
Jr — ieRj'^R. 

Let the scalar function <j)i be in the form 

01 = p'' + p'' = constant, 0'' ~ 

Then omitting higher order terms, equations (16. lip become 
25) aMTj/*; ™ Ti/fc _ n ™ _ „*; 

Let 

mi 



d^^W^, - TUkW^ = 0, mk = p^ak, for fc = 1, 2, 3. 
7712 = Mw, ms = Mz, 
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where Mw and Mz are given by (|6.9|) . Then we derive three vector bosons with 
respective masses as follows: 

which are the same as those derived from the Weinberg-Salam theory. 

Taking divergence on both sides of (I6.1ip and noticing that d'^d^W^^ = 0, we 
obtain that 

(6.16) d^d^(j)\ = mkcjj'l + higher order terms, 

where irik is the constant term in 

a^d^W^ + a^d^A^ — mj, + higher order terms. 

Hence (|6.16p can be regarded as the Higgs field equations, and (l)^ are the Higgs 
particles with mass . 
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